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We perform a systematic search for simple but viable lepton mixing patterns. Our main criterion 
is that the mixing matrix can be parameterized by three rotation angles, which are simple fractions 
of TT. These simple rotation angles possess exact expressions for their sines and cosines, and often 
arise in the flavor symmetry models. All possible parameterizations of the mixing matrix are taken 
into account. In total, twenty successful mixing patterns are found to be consistent with the latest 
neutrino oscillation data (including the recent T2K results) in the CP conserving case, whereas fifteen 
mixing patterns are allowed in the maximal CP violating case. Potential radiative corrections to 
the constant mixing patterns are also calculated by solving the renormalization group equations. 

PACS numbers: 14.60. Pq, 14.60.Lm 



I. INTRODUCTION 

Recent solar, atmospheric, reactor and accelerator neu- 
trino experiments have provided us with compelling ev- 
idence that neutrinos are massive and lepton flavors are 
mixed. In the framework of three-flavor neutrino oscilla- 
tions, the mixing is described by a 3 x 3 unitary matrix 



J 



C12C13 S12C13 

V = I -S12C23 - Ci2S23Sl3e"' C12C23 - Sl2S23Sl3e" 
S12S23 - Ci2C23Sl3e"' -C12S23 



V, which is usually parameterized [l| by three mixing 
angles (6*12, 6'23, and 6*13) and three CP violating phases 
out of which one is the Dirac phase (6) and the other two 
are the Majorana phases (p and a). In the standard pa- 
rameterization advocated by the Particle Data Group 
and in Rcfs. [3[, the lepton mixing matrix reads 



J(5 



si3e 
S23C13 




(1) 



where Sij 



sin 9ij and dj = cos 6ij (for ij = 12,23, 13). 
If neutrinos are Dirac particles, the phases p and a will 
be irrelevant and can be rotated away through a redef- 
inition of the neutrino fields. The latest global analysis 
of current neutrino oscillation data yields [J] 



31.0° < On < 37.1° 

35.7° < 6I23 < 53.1° 

4.1° < 6I13 < 12.9° 



(2) 



at 3(7 C.L., and the best-fit values of three mixing angles 
arc 6I12 = 34.0°, 6I23 = 40.4° and 6*13 = 9.1°. Driven 
in particular by the latest T2K results [5[, ^13 = 0° is 
currently disfavored at the more than 3a level. 

So far it is still unclear how to theoretically understand 
the observed lepton mixing. One tentative way is to start 
with experimental values of leptonic mixing angles and 
conjecture a simple constant mixing pattern, which may 
turn out to be suggestive of the underlying symmetry of 
lepton mixing. In fact, several interesting constant mix- 
ing patterns have been suggested along with the progress 
in neutrino oscillation experiments, and shown to be 



Sl2C23Sl3e C23C13 

I 



derivable from the flavor symmetries. For instance, the 
democratic |6|, bi-maximal [7|, tri-bimaximal [8|, hexag- 
onal 0, [1^ or both golden ratio [il| patterns can be 
realized in models with different discrete flavor symme- 
tries, such as 5*3, A4, S4, A5, dihedral groups, etc., see 
[12| for recent reviews, and [9[ for a summary of proposed 
mixing scenarios. 

Note that all the aforementioned constant mixing pat- 
terns lead to a vanishing ^13, which mainly due to long- 
baseline data is now disfavored at 3cr C.L. Q. While 
more statistics and complementary measurements from 
reactor neutrino experiments will tell us whether ^13 is 
indeed as large as it currently appears to be, it is without 
doubt timely to consider the ways to cope with a sizable 
013 . Indeed, after the results of T2K Q were released, 
several possible ways, with a large range in what regards 
the level of sophistication, to realize a relatively large ^13 
have been discussed in Refs . 113 1. E arlier analyses can be 
found for instance in Refs. |14l. Il5j. 

In this work, we search in a systematic way for suc- 
cessful mixing patterns with initial nonzero ^13. Simple 
requirements are the starting point of our analysis: (i) 



the mixing matrix is a product of three rotations; (ii) 
the angles associated with the rotations are simple frac- 
tions of TT, such that the resulting sines and cosines are 
given by exact expressions. In particular the latter crite- 
rion is reminiscent of the results of many flavor symmetry 
models. We find in total twenty viable mixing patterns 
in the CP conserving case, while fifteen different feasi- 
ble patterns exist in the CP violating case. Furthermore, 
leaving the CP phase arbitrary gives in total 66 successful 
mixing patterns. 

The precision era which neutrino physics has recently 
entered requires that the renormalization effects should 
be taken into account. If the underlying flavor symmetry 
works at some high-energy scale, the mixing angles will 
receive radiative corrections and deviate from the predic- 
tions of a given mixing pattern when running from the 
symmetry scale to the low-energy scale. We will there- 
fore consider the renormalization group equation (RGB) 
effects on our scenarios. 

The remaining part of this work is organized as fol- 
lows. In Sec. ini we recall all different parameterizations 
of a 3 X 3 unitary matrix, and fix our notations. Then, in 
Sec. lIIH we summarize the feasible mixing patterns which 
are compatible with experimental data and in particular 
predict a nonzero 0i3. Scc.|TV]is devoted to a general dis- 
course on the radiative corrections to the constant mixing 
patterns. Finally, we conclude in Sec. El 



II. PARAMETERIZATIONS OF LEPTON 
MIXING MATRIX 

First of all, let us review the classification of possible 
parameterizations of a lepton mixing matrix [l6|. Since 
the Majorana phases can always be recast into a diagonal 
matrix on the right-hand side of V, they have no influence 
on our results and will be ignored for now. If the leptonic 
CP violation is absent, V is simply a 3 x 3 orthogonal 
matrix and can be written as a product of three rotation 
matrices with three different rotation angles ("iJi, i?2, ^3), 
i.e. 

V = R,,idi)Rkli^2)Rmni^3) , (3) 

where ij, kl, mn = 12, 23, 13 and 



cos-i? sin-i? 0^ 

— sin -d cos "d 

Of, 

f 

cost? sint? j , (4) 

— sin d cos d j 

cost? sini?^ 
I 

- sin 1? cos 1? I 



Note that the order of the rotations is not specified. The 
Dirac-type CP violating phase ip can be included in the 




above parameterization by replacing the entry "1" with 
a phase factor e~^'^ in the second rotation matrix on the 
right-hand side of Eq. dH). Taking Ri2{'d) in Eq. (4) for 
example, we have 



Ri2{^,v>) 



cos 1? sin -d 
— sin d cos ^ 



(5) 











Although there are several different ways to introduce the 
CP violating phase, the choice in Eq. ([5]) is advantageous 
in the sense that the phase parameter ip is always located 
in a 2 X 2 submatrix of V, in which each element is a sum 
of two terms with the relative phase ip. 

In Ref. [la] it was shown that only nine distinct pa- 
rameterizations exist, namely 



Pi 

P2 

Ps 
Pa 
P5 
Pe 
Pi 
Ps 
P9 



V 
V 
V 
V 
V 
V 
V 
V 
V 



Rl2i'&l)R23i'&2 
R23{^l)Rl2{^2 
R23{^l)Rl3{^2 
Rl2{^l)Rl3{^2 
Rl3{^l)Rl2{^2 
Rl2{^l)R23{^2 
R23{^l)Rl2{^2 
Rl3{^l)Rl2{^2 
Rl3{^l)R23{^2 



V)Rl2{^^) . 
^)R23\^3) , 
^)Rl2i^3) , 
^)R23\^3) > 
V^)Rui^3) , 
^)Rl3i^3) , 
^)Rl3(^3) , 
^)R23i^3) , 
^)Rl2i^3) ■ 



(6) 



Here i?,"-^(i?) 



Ri 



^ \"/ -Hj(— I?)- Three of the nine parameteri- 
zations belong to the class ij = mn ^ kl and six to the 
class ij ^ kl ^ mn. Note that P3 is just the standard 
("PDG") parameterization in Eq. ([IJ, up to a simple 
phase redefinition. 

The effects of CP violation are usually characterized 
by the Jarlskog invariant Jqp |f 7| , which is defined as 



Im[K. 



y^j^aj^iSi 



]=JcP E 



E 



{£ijk£al3-t) ■ (7) 



7=e,/i,r fe=l,2,3 



It is straightforward to verify that the Jarlskog invariant 
is given by 



-^CP = SiCiS2C2S3C3Smip 



for Pi , P2 and P5 , and 



-^CP = siCiS2C2S3C3sin(p 



(8) 



(9) 



for P3,P4,Pq,Py,Ps and Pg, where Si = siwdi and q = 
cos-i^i for i = 1,2,3. 

It is worthwhile to remark that although these nine 
parameterizations are mathematically equivalent, one of 
them may turn out to be more useful than the others for 
a specific problem. For instance, three mixing angles in 
the standard parameterization P3 can be unambiguously 
extracted from neutrino oscillation experiments, which 
is not the case for the other parameterizations. The fla- 
vor symmetry behind the observed lepton mixing pattern 
may be manifest in a certain parameterization, which is 
currently unknown to us, so we consider all possibilities 
in Eq. ^ when searching for viable constant lepton mix- 
ing patterns. 
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5 + V5 
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2 + ^3 
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TABLE I: Exact expressions of sin^ i? and cos^ -ff for simple 
rotation angles i? = vr/n with n = 1, 2, 3, 4, 5, 6, 8, 10, 12. 



the standard parameterization defined in Eq. ([T]) and 
Pj(ni, 712, na), one can immediately extract the three 
"standard" leptonic mixing angles 



sm 



712 = 



sm 6/23 



sin 6*13 



^e2\ 



"631 

I-IK3I 
^e3|^ 



IK: '2 



(12) 



III. SIMPLE ROTATION ANGLES AND 
VIABLE MIXING PATTERNS 

One immediate question arises: which values of the ro- 
tation angles should we use to account for the observed 
lepton mixing? To make the mixing patterns simple and 
suggestive of flavor symmetries, we set the following cri- 
teria: (i) three rotation angles are simple fractions of tt, 
i.e. ir/n with n being an integer; (ii) the choice of n is 
governed by the requirement that the sines and cosines of 
the rotation angles possess exact expressions, i.e. they are 
expressible as simple terms involving only square roots. 
Although these two criteria are essentially set to make the 
lepton mixing matrix as simple as possible, they might 
be realized in the flavor symmetry models (e.g., the di- 
hedral group Dn [18| ) , where the lepton mixing patterns 
can be intimately related to Clebsch-Gordan coefficients 
favoring small integers and their square roots. 

Up to n = 12. one has the following nine rotation an- 
gles satisfying the above criteria 



-& 



{TT TT TT TT TT TT TT TT ^ 
^'2' 3' 4' 5' 6' 8' 10' 12/ 



(10) 



for which the values of sin^ 1} and cos^ 'd have been listed 
in Table m Note that for n > 12, there also exist simple 
angles whose sines and cosines have exact expressions. In 
general, one can always implement the relation 



cos — 
2 



1 — cos d 



(11) 



to calculate cos(7r/2n) if cos(7r/ri) is already known. 
However, it should be rather difficult to relate the expres- 
sions with three or more square roots to a flavor symme- 
try, and in addition the resultant lepton mixing pattern 
becomes very complicated. Furthermore, as n increases, 
"d becomes smaller, and because the smallest mixing ma- 
trix element Ves is generally related to the smallest rota- 
tion angle, the resulting mixing matrix ends up not being 
in agreement with experimental data. 

The lepton mixing matrix V can be obtained by in- 
serting the rotation angles di = TT/rii (for z = 1, 2, 3) into 
Eq. ©, where n, e {1,2,3,4,5,6,8,10,12}. For later 
convenience, we denote the mixing pattern constructed 
through the parameterization Pj (for j = 1,2,..., 9) 
with rotation angles di = Tr/n^ as V = P,(ni, 712,713) 
for a given CP violating phase ip. Comparing between 



which allow us to confront the lepton mixing matrix 
Pj(7ii, 7i2, 713) with neutrino oscillation data and thus find 
out the viable patterns. Note that with our set of rota- 
tion angles we can not generate a viable pattern with 
the standard parameterization P3 in Eq. ([1]), because 
jVesI = sin6'i3 in this case and the smallest possible an- 
gle in our scenario is 7r/12, implying a too large value of 
sin^ 6113 = 0.07. The largest allowed value of 6*13 approx- 
imates to 7r/14, which does not possess an exact expres- 
sion for its sine or cosine. 

There exist nine distinct parameterization schemes and 
9x9x9 different combinations of rotation angles in each 
scheme, so we are left with 9"^ = 6561 possible mixing ma- 
trices for a fixed CP violating phase ip. We have numeri- 
cally studied all these possibilities for the CP conserving 
case if = and the maximal CP violating case ip = tt/2 
by comparing the predicted mixing angles with the 3ct 
global- fit data given in Eq. ([2]). The viable cases arc listed 
in Table HT] We find that 20 mixing patterns are allowed 
in the CP conserving case, while only 15 patterns are 
compatible with experimental data in the maximal CP 
violating case. All in all, the feasible patterns for differ- 
ent CP violating phases are different. Only ^7(5, 5, 12) 
is allowed for both (p = and ip = tt/2, showing the 
importance of the CP violating phase in searching for 
successful lepton mixing patterns. We should note here 
an important difference between the standard parame- 
terization P3 in Eq. ([1| and the others. In the standard 
parameterization, a choice of three angles can directly be 
confronted with experiments, namely knowing | I42I, \Ve3\ 
and I V^3 1 fixes the mixing angles and leaves the phase un- 
determined. Consider now P2, whose explicit form is 



V 



Cl 


S2C3 


-S2S3 


-CIS2 


C1C2C3 + sisse""^ 


-C1C2S3 + S1C3' 


SIS2 


-S1C2C3 -|-ciS3e""^ 


S1C2S3 -l-cicse 



-\ip 



Since |Ve2|, jV'esI and \V^3\ are known experimentally, two 
angles (^2 and 2^3), as well as some combination oi-di and 
the CP phase ip, can be determined. On the other hand, 
choosing three angles, as we have done in this paper, fixes 
the CP phase ip at the same time. This is true for all the 
alternative parameterizations Pi, 2. 4, 5. 6, 7. 8, 9 of the lepton 
mixing matrix. 

We continue by recommending a few simple but inter- 
esting mixing patterns and discuss their implications for 
the leptonic mixing angles and leptonic CP violation: 



Patterns 


V3 = 


V3 = 7r/2 1 


612 


02-A 


^13 


6*12 


623 


0l3 


Jcp i%J 


Pi(10,4,4) 


32.1° 


43.6° 


12.6° 










Pi(10,4,5) 








36.9° 


43.6° 


12.6° 


4.9 


Pi(10,4,6) 








31.6" 


43.6" 


12.6° 


4.5 


Pi(12,4,4) 


34.3" 


44.0° 


10.5° 










Pi(12,4,5) 








36.6° 


44.0° 


10.5° 


4.2 


Pi(12,4,6) 








31.1° 


44.0° 


10.5° 


3.8 


P2(3,5,10) 


34.6" 


45.3° 


10.5° 










P2(3,5,12) 


35.1" 


47.8" 


8.8" 










P2(4,5,10) 








34.6° 


45.0° 


10.5° 


4.1 


P2(4,5,12) 








35.1" 


45.0" 


8.8" 


3.5 


P2(5,5,10) 








34.6" 


37.2" 


10.5" 


3.9 


P2(5,5,12) 








35.1° 


36.8° 


8.8° 


3.3 


^4(6,8,4) 


36.3° 


48.9° 


6.8° 










P4(6, 10,4) 


33.4° 


47.0° 


9.5° 










P4(6,10,5) 


34.3° 


39.5° 


4.4° 










P4(6,12,4) 


31.5" 


45.9" 


11.2° 










P4(6,12,5) 


32.7° 


38.1° 


6.5° 










P4(8,6,5) 


36.0° 


44.9° 


8.6° 










P4(8,6,6) 


35.1° 


39.9° 


12.0° 










P4(10,6,4) 


34.0" 


52.0° 


6.8" 










P4(10,6,5) 


32.7" 


44.3° 


11.7° 










A(12,6,4) 


32.1" 


51.7" 


9.1° 










P5(5,4,3) 


35.7° 


38.7° 


11.6° 










P7(4, 5, 12) 








36.9" 


45.0" 


12.1° 


4.8 


Pr(5,5,12) 


36.9° 


45.0° 


12.1° 


36.9° 


36.4° 


12.1° 


4.6 


P7(6,5,12) 


36.9° 


39.0° 


12.1° 










P9(10,4,5) 








36.9° 


46.4° 


12.6° 


4.9 


P9(10,4,6) 








31.6° 


46.4° 


12.6° 


4.5 


P9(10,4,8) 


35.4" 


46.4° 


12.6° 










A(12,4,5) 








36.6" 


46.0" 


10.5" 


4.2 


A(12,4,6) 








31.1° 


46.0° 


10.5° 


3.8 


P9(12,4,8) 


33.2" 


46.0" 


10.5" 










P9(12,5,5) 








36.4° 


36.9° 


12.1° 


4.6 


P9(12,5,8) 


31.5° 


36.9° 


12.1° 











TABLE II; List of viable mixing patterns, in which the lep- 
tonic mixing angles are within the 3(t ranges of the global-fit 
data. In the case oi tp = 7r/2, the Jarlskog invariant Jcp is 
given in the last column. 

(1) Pattern Pi (12, 4, 4) with i^ == - The lepton mixing 
matrix takes the form 



(13) 



/V3-1 1 v^+l 
4^2 ' 4 

V — v^+1 ^-1 

4^2 4 

V ~5 


v^-l %/3+l 
4^2 4 

v^+l 1 V^-1 
4^2 ' 4 

~2 


\/3-l 

4 

\/3+l 

4 

1 

%/2 


which leads to 






sin^^ia = 


-li^-^] 


; 


sin^023 = 


2 + V3 
6 + 73 ' 




sin^ 9i2 = 


1 V2 


; 



(14) 



6 + ^/3 

or explicitly 6*13 « 10.5°, 6*23 « 44.0° and 6*12 « 34.3°. 
Such a mixing pattern is in excellent agreement with neu- 



trino oscillation data. The predictions for 612 and ^23 
fall into the la ranges, while that for 613 is even slightly 
larger than the best-fit value but well within the 2a range. 
(2) Pattern P4(12, 6, 4) with ip = - The lepton mixing 
matrix takes the form 



v = 


/ 3-I-V3 
4%/2 
3-\/3 
4\/2 


3\/3-l 3-%/3 

8 8 

3+V3 3V3-I-1 

8 8 

V3 v^ 

2^2 2v^ 


which leads to 




sin^ 6*13 = 


^ l(-v^) 


sin2 6'23 = 


14 + 3a/3 
26 + 3\/3 ' 


sin^ 012 = 


14-3V3 
26 + 3\/3 ' 


or explicitly 6*13 « 9.1°, 


6*23 ~ 51.7° ai 



\ 



(15) 



(16) 



12 « 32.1°. 

Note that the prediction for 6*13 in this mixing pattern is 
almost the best-fit value. 

(3) Pattern Pi (12, 4, 6) with ip = tt/2 - The leptonic 
mixing matrix takes the form 



v = 


/v/3-l-i%/2(3+y3) 

Q 

\/3+l+i%/2(3-%/3) 
8 

V 2V2 


3-\/3-l-iV2(V3+l) 

Q 

3+v^-i\/2(\/3-l) 

2V2 


\/3-l 
4 

Vs+i 

4 

1 

V2 


which leads to 






sin^0i3 = 


-li^-^]^ 




sin^6'23 = 


2 + Vi 
6 + V3 ' 






sin^eiis == 


10 -V3 





,(17) 



(18) 



24 + 4^3 ' 



or explicitly 6*13 « 10.5°, 6I23 ~ 44.0° and 6I12 w 31.1°. In 
addition, the Jarlskog invariant is Jqp = •\/6/64 « 3.8%. 
(4) Pattern P9(12,4,6) with ^ == 7r/2 - The lepton 
mixing matrix takes the form 

-y3+l-i\/2(3+\/3) -3+\/3-l-i\/2(V3-H) Vlz-l^ 



(19) 



8 
\ -y3-l+i\/2(3-\/3) 


2\/2 
-3-\/3-i%/2(\/3-l) 


4 

1 

V3+1 


^ 8 


8 


4 


which leads to 






sin2 6'i3 = 


li'-^)' 




sin^6'23 = 


4 




6 + V3 ' 




sin2 6'i2 = 


10 -V3 





(20) 



24 + 4\/3 ' 



or explicitly ^ig w 10.5°, 6*23 « 46.0° and 6*12 w 31.1°. In 
addition, the Jarlskog invariant is Jqp = •\/6/64 sa 3.8%, 
which is the same as that in the previous case. Note that 
the mixing patterns in Eqs. (17) and (19) differ only in 
the predictions of 623, i.e. 6*23 = 44.0° in the former case 
while 6*23 = 46.0° in the latter. 

As mentioned before, the CP violating phase plays an 
important role in searching for successful lepton mixing 
patterns. In this regard, we have so far focused on the 
special cases with (p = and (p = tt/2. Choosing an ar- 
bitrary CP violating phase, i.e. (p G [0,tt], leads to the 
66 viable patterns in Fig. [TJ One can immediately con- 
struct the lepton mixing matrix from the Pj(ni, 712, 713) 
notation. It is amazing that from the 6561 possible mix- 
ing patterns, only about 1% are compatible with current 
oscillation data. 



Ve^yl^Vr M), 



IV. RADIATIVE CORRECTIONS 



Now we proceed to consider possible radiative correc- 
tions to the mixing patterns in Table |lll ^ As argued in 
SecU the mixing patterns under consideration may arise 
from certain flavor symmetries preserved at high-energy 
scales, such as the grand unification (e.g. A ^ 10^® GeV) 
or the seesaw scale (e.g. A ^ lO^'' GeV), whereas the lep- 
tonic mixing parameters are determined or constrained in 
neutrino oscillation experiments at low energies. The gap 
between the high-energy predictions and the low-energy 
measurements is bridged by the RG evolution, which may 
significantly change the model predictions. On the other 
hand, the RG running effects could also serve as an expla- 
nation for the discrepancy between the fiavor symmetric 
mixing pattern and the observed one. 

The RGBs for leptonic mixing parameters have been 
derived within various theoretical frameworks [l9[. In 
the supersymmetric theories with large tan/3, it has been 
found that the RG evolution may lead to significant mod- 
ifications to the mixing parameters, in particular the so- 
lar mixing angle 6112 (see e.g. Ref. [20| and references 
therein). To be explicit, we write down the RGEs for 
three leptonic mixing angles in the approximation of r- 
lepton dominance (i.e., Y^Y^ sa diag{0,0, y^} in view of 



^ In the fiavor symmetry models, the corrections to fermion masses 
and mixing patterns may also originate from the flavor symmetry 
breaking and the higher-dimensional operators. Since the signif- 
icance of these potential corrections is highly model-dependent, 
we shall concentrate on the generic RGE corrections in the cur- 
rent study. 



' g^^Am^ ^™i + ""2 + 2miTO2C2(p-.)J , 

^?^T^12'^12^23'"23'^3 



- (1 + C) m2C(^2a+S) - Cm^cs] , (21) 

'23 - -V^hr^ ^ H2 {ml +ml + 2m2m3C2a) 



where 9ij — d9ij/dt with t = ln(/i//io), C = 
AmL/Am^jj^ with Am^^, = m'^ - m^ ^ 7.6 x 10"^ eV^ 



'sol 



sol 



and \Aml^^\ = |m§ - to^| w 2.3 x 10"^ eV^ at the low- 
energy scale, and yr denotes the Yukawa coupling of tau 
lepton. In the standard model (SM) C ~ —3/2 while 
in the minimal supersymmetric standard model (MSSM) 
C = 1. Terms of 0(^13) have been safely neglected in 
Eq. ([2T|) . where we have defined C2(p-a) = cos2(p — a), 
C2p = cos 2p and so on. 

As obvious from the above beta functions of the mix- 
ing angles, 9i2 receives typically larger RG corrections 
than 6^23 and ^13 , whose corrections are of the same or- 
der. Furthermore, when running from a high-energy scale 
to the electroweak scale Aew = 10^ GeV, the radia- 
tive corrections to 9i2 are typically (see below) positive 
in the MSSM, i.e. 6*12 (A) < 6*12 (Aew)- In contrast, in 
the SM 6112 receives only negative corrections because of 
the sign flip in C. However, the RG effects in the SM 
are generally small due to the absence of tan /3 enhance- 
ment. As for 023, the RG corrections could be either 
positive or negative, depending mainly on the model and 
the neutrino mass ordering. In the normal mass order- 
ing (toi < 7722 < "^3) both typically decrease in the SM 
and increase in the MSSM, whereas for the inverted mass 
ordering (7713 < mi < 777,2) the behavior is opposite. Fi- 
nally, nonzero ^13 can run in both directions, and receives 
corrections of the same order as ^23 • 

When searching for viable lepton mixing patterns, we 
have ignored the Majorana phases, which indeed do not 
change the mixing angles. However, they may play a 
significant role in the RG evolution of mixing angles, 
in particular for 9i2. In the CP conserving limit with 
p = cr = 0, 012 c>c (7771 4-7772)^/A77igQ[; which may strongly 
boost the RG running in the nearly degenerate or in- 
verted mass ordering. On the contrary, in the limit of 
p — a = tt/2, 9i2 (X {1712 — 77ii)/(77i2 + TUi), wliicli is 
always smaller than one and hence suppresses the RG ef- 
fects on 012. In the latter case, the running mainly comes 
from the next-to- leading-order terms of 0i3 , and thus 0i2 
may run to a slightly smaller value. Since both p and a 
are entirely unconstrained in oscillation experiments, we 
shall allow them to freely vary between and n. 

We have numerically solved the full set of RGEs for 
leptonic mixing angles. More explicitly, the lepton mix- 
ing patterns in Table [ll] arc assumed at a cutoff scale 
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FIG. 1: Viable patterns for an arbitrary CP violating phase ip. The red, green, and blue bars denote 
the allowed ranges of ip, for which the predicted mixing angles fall into their la, 2a, and 3a intervals, 
respectively. Patterns with ip — are written in blue, while patterns with ip = 7r/2 are written in red. The 
pattern accommodating both ip = and ip = n/2 is highlighted in green. 
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FIG. 2: RG correction to 612 in tlie SM (green bars), and 
in the MSSM (red bars for the normal mass ordering and 
blue bars for the inverted mass ordering). The initial values 
are labeled as black diamonds. Note that <^ = is assumed 
and tan/3 = 10 is adopted in the MSSM. Furthermore, we 
allow the lightest neutrino mass mi (or ma) to vary in the 
range of (0...0.1) eV in the MSSM, whereas mi (or ms) 
varies in the range of (0 . . . 0.2) eV in the SM. The vertical 
lines correspond to the best-fit value and the Icr, 2a and 3cr 
intervals. The mixing scenarios written in red are not valid 
without RG corrections. 



A = 10^° GeV, and then we evolve the mixing parame- 
ters down to the electroweak scale Aew in order to com- 
pare them with experimental data. Note that the depen- 
dence of RG corrections on the cutoff scale is logarithmic 
and therefore the precise value of A is not quite rele- 
vant. Other physical parameters, e.g. gau ge couplings 
and fermion masses, arc taken from Rcf. [23] • In the case 
with (p = 0, we show in Fig. [2] the RG evolution of 6*12 
in the MSSM with tan/3 = 10 as well as in the SM. One 
can observe that in the MSSM some viable mixing pat- 
terns can in principle receive RG corrections so large that 
they are no longer valid, for instance ^9(10, 4, 8) in the 
inverted mass ordering. However, as mentioned above, 
by choosing p~ a = 7r/2 one can suppress the running of 
012, a situation with interesting consequences for neutri- 
noless double beta decay, as discussed in Ref. [l4|. Note 



FIG. 3: The RG corrections to 613 in the SM (green bars), 
and in the MSSM (red bars for the normal mass ordering and 
blue bars for the inverted mass ordering). The initial values 
are labeled as black diamonds. The input parameters are the 
same as in Fig. [2] Note that the corrections in the SM appear 
to be invisible. The mixing scenarios written in red are not 
valid without RG corrections. 



that, although the predictions on 9i2 arc same for both 
P2(3, 5, 12) and P4(8, 6, 6), they suffer from different RG 
corrections, reflecting the importance of 6*23 in the RG 
evolution [cf. Eq. ((2T|) ]. In the SM, the RG corrections 
are in general small, and all the viable mixing patterns 
remain valid at the electroweak scale. In Fig. [5] we have 
also shown twelve new mixing scenarios, which generate 
too small or large ^12 when not corrected by RG effects, 
but can enter the allowed range after RG corrections. 
These cases are indicated by writing them in red (see 
below). 

For completeness, wc also show the RG corrections to 
^13 in Fig. |3l As expected, no visible effects can be seen 
in the SM, while in the MSSM, the RG running may 
lead to tiny deviations (less than one degree) from their 
initial values. Even in this case, the values of ^13 from 
some mixing patterns listed in Table |ll] may exceed their 
3(7 ranges, e.g., the mixing pattern Pi(10,4,4). Further- 
more, the RG effects on ^23 are of similar size to those on 
6*23 [cf. Eq. (pij) ]. and therefore no significant radiative 
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FIG. 4: The RG corrections to ^12 in the SM (green bars), 
and in the MSSM (red bars for the normal mass ordering and 
blue bars for the inverted mass ordering). The initial values 
are labeled as black diamonds. We use the same parameters 
as in Fig. [2] except <p = 7r/2. 



corrections can be acquired. 

Now we consider the maximal CP violating case with 
ip ~ 7r/2. The RG corrections to Q\i are depicted in 
Fig. m Similar to the CP conserving case, sizable RG 
effects can be present. In the MSSM and in the case of the 
inverted mass ordering, one can observe that for instance 
Pg(12,5,5) can become incompatible with experimental 
data, and similarly for Pg(12,4, 6) in the SM. On the 
other hand, there are eleven new mixing schemes which 
become valid only after sizable RG corrections. 

Let us turn to the mixing patterns which are not com- 
patible with data at the cutoff scale but can be modified 
to the proper parameter ranges with the help of the RG 
running. This is in particular possible for the mixing 
patterns with both ^23 and dn in the currently-favored 
ranges but a smaller Q\2^ since dyi may be lifted up into 
the correct parameter interval in the MSSM. For exam- 
ple, the mixing patterns with Qyi = 30°, or sin0i2 = 1/2, 
could be of particular interest [a,[l3|^- We have also per- 
formed a systematic search for such mixing patterns in 
our scenario. As mentioned above, in the case of (/? = 



((/9 = 7r/2), there exist twelve (eleven) new patterns. Now 
we shortly discuss three of them: 

(1) Pattern ^5(4,4, 3) with (p = - The lepton mixing 
matrix is 



(22) 
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/v^+i 
k 


1 V2-V3 
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sin^ 6*13 — 


^(5-2^/6) 
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sin^6'i2 = 


II + 2V6 ' 
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11 + 2V6 ' 



(23) 



or explicitly 6*13 « 4.6°, 6*23 « 38.0° and 6*12 ~ 30.1°. 
Evolving the mixing angles through the full set of RGEs 
in the MSSM with tan/3 = 10, we obtain, at the elec- 
troweak scale, 6*12 G [30.1°, 36.0°] for the normal mass 
ordering and 6*12 S [30.1°, 36.9°] for the inverted mass 
ordering [cf. Fig. [2]. Thus the mixing angles become 
consistent with the experimental data with the help of 
radiative corrections. 

(2) Pattern ^2(3, 6, 12) with ip = - The lepton mixing 
matrix takes the form 



v = 


/ VS \/3+l v^-1 
2 4%/2 4%/2 
1 9-V3 3(%/2+%/6) 
4 8%/2 16 
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14 + \/3 ' 



(25) 



or explicitly 6*13 « 7.4°, 6*23 « 46.9° and 6*12 « 29.1°. 
Again we evolve the mixing angles via the RGEs in the 
MSSM with tan ^ = 10, and obtain 6I12 £ [29.2°, 37.7°] in 
the normal mass ordering case while 612 € [29.1°, 38.8°] 
in the inverted mass ordering [cf. Fig. [2]. These values of 
6']^2 are well compatible with the 3cr range in Eq. (2). 

(3) Pattern ^2(4,6, 12) with (p ^ tt/2 - The lepton 
mixing matrix takes the form 



^ Another example of this kind is the so-called tetra-maxinial mix- 
ing pattern with 6*12 f» 30.4° (tan 612 = 2 - V2) [H. 
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which leads to 
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T^(2-V3) 
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2 + V3 



(27) 



14 + V3 ' 



or explicitly 6I13 w 10.5°, 623 == 45.0° and 6*12 ~ 29.1°. In 
addition, the Jarlskog invariant is Jcp = •\/3/64 ps 2.7%. 
Similar to the above two patterns, the RG running in 
the MSSM leads to 6*12 G [29.0°, 36.6°] for the normal 
mass ordering and 9i2 £ [29.0°, 37.1°] for the inverted 
mass ordering [cf. Fig. |4], which are in agreement with 
the experimental data. 

In general, the evaluation of RGEs may play a crucial 
role in searching for realistic mixing patterns. Especially 
in the supersymmetric case, a larger tan/3 typically leads 
to more significant RG corrections. In this sense, the 
mixing patterns with smaller ^12 could be more favor- 
able for a larger tan/3. It is also worthwhile to stress 
that the RGEs under discussion are given in the effec- 
tive theory approach, which is essentially the same for 
different kinds of flavor or seesaw models. In the realistic 
flavor symmetry models, the flavons might induce addi- 
tional contributions to the RGEs. If a seesaw model is 
considered, the RG running between the seesaw thresh- 
olds may also lead to remarkable modifications, and thus 
should be carefully treated [2J|. A thorough survey on 
the RGEs in a specific model and on the different values 
of tan /3 is beyond the scope of this work, and we refer 
the readers to Refs. [l9l-[2ll for more detailed discussions. 



V. CONCLUSIONS 

Motivated by the recent indications of a nonzero ^13, 
we have performed a systematic search for simple but vi- 



able lepton mixing patterns by setting two criteria: (i) 
the lepton mixing matrix is parameterized by three rota- 
tion angles, which are simple fractions of tt; (ii) the sines 
and cosines of these rotation angles possess exact expres- 
sions. In total, we have found 20 viable mixing patterns 
in the CP conserving limit, while 15 viable patterns exist 
in case of maximal CP violation. Moreover, in the most 
general cases with the CP phase unconstrained, only 66 
mixing patterns out of 6561 combinations are found to 
be compatible with current data. 

Furthermore, radiative corrections to the mixing pat- 
terns have been calculated by solving the RGEs of lep- 
tonic mixing parameters. We have shown that the RG 
running can induce sizable corrections to the lepton mix- 
ing patterns, which eventually could render some pat- 
terns to be unsuccessful in describing lepton mixing. On 
the other hand, we have also pointed out some interest- 
ing mixing patterns, which are incompatible with current 
oscillation data at the high-energy scale but become vi- 
able at the low-energy scale after the RG corrections are 
properly taken into account. 

We hope that the successful constant mixing patterns 
found in this work can be helpful in searching for the 
underlying flavor symmetries and shed some light on the 
final solution to the flavor puzzle. At least, they could 
serve as a useful phenomenological description of lepton 
mixing. 
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